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Abstract 

The following questions, concerning to the application of the harmonic oscillator represen- 
tation (HOR) in the theory of scattering and reactions, are discussed: the formulation of the 
scattering theory in HOR; exact solutions of the free motion Schroedinger equation in HOR; 
separable expansion of the short, range potentials and the calculation of the phase 1 shifts: 
’’isolated states” as generalization of the Wigner-von Neumann hound states embedded in 
continuum; a nuclear coupled channel problem in HOR; the description of true three body 
scattering in HOR. As an illustration the soft dipole mode in the 11 Li nucleus is considered 
in a frame of the 9 Li + n + n cluster model with taking into account of three body continuum 
effects. 


1 Introduction 

Usually harmonic oscillator wave functions are used for the description of bound states of quantum 
systems that belong to the discrete spectrum [1]. In this talk the application of the harmonic 
oscillator (HO) basis to the solution of the scattering problem, i.e. in continuum , will be discussed. 

This line of investigations was begun in Refs. [2] (see also the papers cited there) and inde- 
pendently in the papers of Kiev [3] and Moscow [4] groups. The similar approach, also connected 
with an application of the HO basis to the scattering problem, was developed by the Hungarian 
group [5], 

In order to illustrate the essence of the approach to the scattering problem in the harmonic 
oscillator basis, we shall consider at first the simplest problem of the scattering of a single particle 
by the central potential U(r) [4]. Thus, we come to the Schroedinger equation 

+ V{r) \ V’im(r) = £ tl’im ( r ) . ( I ) 

Its solution V’bn(r) = £ti(r)Yt m {Q) will besought in the form of an expansion in the eigenfunctions 
of the harmonic oscillator 

fl/(r) = t^niRniir) (2) 

71 — 0 
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r /£(+!/2( r 2 ) e -r 2 /2 


( 3 ) 



where 


Re(r) = (-l) n 


2n! 


r(n + / + 


§)) 


1/2 


is the radial wave function of three-dimensional harmonic oscillator. This wave function cor 
responds to the eigenvalues of the oscillator energy E™ c = (2n + / + 3/2)hu> [1]. The value 
r Q = (h/muyl 2 is selected as a length scale in relations (1) and (2). Here u is the oscillator 
frequency; the energy e = q 2 / 2 is measured in units hu>] the wave vector k is expressed in units 
r -\ q = kr 0 is the dimensionless momentum. Substituting expansion (2) in (1) and multiplying 
(1) scalarly by Rni^imi we obtain the following equation determining the wave function i/’jm the 
harmonic oscillator representation (HOR): 


£(#w - 6 nn ,e)C n .i = 0 , n = 0,1,2,.... (4) 

n * 


Here, H = T + V and only the following matrix elements of the kinetic energy operator T 
are nonvanishing: 


1 nn — 1 


l_ 

2 


T = 

s nn 


2 


n(n + (+ j) 

( 2 » + l + j) , 


) 


Tnn + 1 


\_ 

2 



n + l + 



-.1/2 


p 2 ! 2 


( 5 ) 


As to behaviour of the coefficient C n t for n > N , their asymptotics are similar to the asymptotic 
of the wave function in the coordinate representation [3] if r is substituted by 2n 1 / 2 r 0 : 


C ni ~ 2n 1/4 rp lm (2y/nr 0 ) , n -> oo . (6) 

this result can be obtained if the WKB expression for the oscillator function is substituted 

in the expression for the coefficients 


Cnl =< V’nlm(r)|t/7m(r) > (7) 

and the integral (7) is calculated by the stationary phase method. The result (6) follows also from 
the fact that the finite-difference equation 


n n 


+ / + 


5) 


1/2 


C n -u + ^2n + / + - — q 2 ^j 


C nl 


(n + 1) (n + / + - 


- 1/2 

Cn+11 T 2 < n/|V|n / > C n <i = 0 

n' 


( 8 ) 
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in the limit n >■ v = 1/2 + 3/4 can be replaced by the following second-order differential equation 
[4]: 

X[' — — ~t ^ ■-A’/ ~ / R(x, x')\/xx'Xi{x')dx + qf 2 X/ = 0 . (9) 

X 2 io 

Here x = 2(n + i/) 1/2 , X/(x) = x xl2 C nl . The Eq. (9) should be solved at the boundary condition 

Xi(2y/^J) = 0 . ( 10 ) 

Thus, in the asymptotic limit for large n, the wave function of our system Xi for the partal 
wave with angular momentum / in the HOR obeys the conventional Schroedinger equation with 
nonlocal potential 


V(x,x')\fxx' ~ 2 < nl\V\n'l > [(n + i ')(«' + v)] 1 * 4 (11) 

where the value 2(n + i/)^ 2 ro plays the role of "coordinate”. In actual calculations, the potential 
matrix has to be cut off by the condition V nn < = 0, if n or (and) n' > N. Then a set of equations 
(4) can be splitted in two parts: 


N 

a) n < N , r. ( H nn ' — £$ nn ')C n >i = — <5 n ^7yvAf+iC’w + i/ , (12a) 


n '=0 

b) n > N + \ , T nn -\C n -u + ( T nn — e)C n i + T nn+1 C n+ u = 0 . ( 1^^) 

Thus, the coefficients C nl with n > N obey the equation of free motion (12b) or, in the asymptotic 
limit of continuous n, the Schroedinger equation of free motion 


X 


ft 

l 


1(1 + 1 ) 


Xi + q 2 X, = 0 . 


It means that the condition 

C n/ ~ 2n 1/4 e- 2 ^ (13a) 

(where e = —it 2 / 2 is the binding energy) must be satitfied for the bound states. The coefficients 
C nt for the scattering problem have the following asymptotic behaviour: 

C n i ~ 2n 1 ^ 4 s\n(2q\/n — /7r/2 + Si) ( 1 36) 


where e - q 2 j 2. According to Eq.(6) the phase shift 6, in Eq. (13b) coincides with the standard 
phase shift of the wave function in coordinate spase. For the decaying resonance states, we get 
(see in [3]): 

C„, ~ 2n^e^ . (13c) 


If the calculations are made up to sufficiently high values of Af > 1 it is possible to use the 
asymptotic expressions (13) [3] . At modest N it is necessary to use the exact, rather than 
approximate, solution for the equation of free motion (12b) which was found in Refs. [2, 4] in 
order to calculate the binding energy, the scattering phases etc. Before considering the solution 
for the equation of free motion, we shall note that the solution for the set (12) is equivalent to the 
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solution for the Schroedinger equation with Hamiltonian H = T + V N containing the many-term 

separable potential 

N 

V N = ^2 < nl\V\n'l > \nl X n'l\ (14) 

n y n f 

with harmonic oscillator formactors. The technique of solving such an equation in the frame of the 
momentum representation was described in [5, 6]. Here we shall describe an alternative method 
for solving the same problem in HOR. 


2 Solution for the equation of free motion in the har- 
monic oscillator representation 

Consider first the case of positive energy £ = q 2 / 2 > 0. The Schroedinger equation of free motion 
in the coordinate space has two linear-independent solutions (regular and irregular) [4]: 


- *(*»•) 


ft 


irreg 


ni(kr) 


1 • n 

— S ,n( fc - j) , 

1 lb /7r \ 

' - — cos(kr ) . 

kr 2 


5) 


In accordance with this, the finite difference equation of free motion (12b) will have also two 
fundamental solutions in the HOR [4] namely the regular solution 


C%’(q) = 


2 r(» + / + t ) 


1/2 


r(n + i) J r(/ + 1) 




-n,l+--q 2 ) = 


03/2 

z 


= (-!)”/?„/(</) ~ —^n l/4 j,(2y/nq) 

satisfying the boundary condition (6) CZu = 0, and the irregular solution 


16) 




' 2i> + i) 
r(n + / + §) 


1/2 


(-!)'</ 




,-? 2 /2 


r(-/ + i) 

23/2 

1 / 4 , 


M( — n — / — </ 2 ) — 


7T 


1/2 


“71 7 71/ 


/(2v/n?) 


.(17) 


which is singular at the point n = — 1. 

The Casorati determinant K n i for these two solutions which plays the same role for the differ- 
ence equations as the Wronskian for the differential equations [7] is of the form: 


I\ n i = T y 


n + ln 


s^reg x^irreg 
ill ^~'nl 

fittg ptrreg 

L 'n+1/ °n+l/ 


nq 


(18) 


Since K n i 7^ 0 for any values of n and /, the expressions (16) and (17) constitute the fundamental 
set of solutions for equation (12b). A11 arbitrary solution for (12b) may be presented as a linear 
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combination of fundamental solutions. 


In particular, the solution for the set (12) for n > N must 


be of the form 

c nl ( q ) = co S s l c:: l 3 ( q )+ s \ns l c:;r 9 (q) 


(19) 


whence it follows that 


tan Si = 


s t s-treg ^ r^rtg 


s' i srirreg s~i syirreg 
^ nA ' n +\l — ^ n + U^ n ( 


( 20 ) 


The equivalent pair of fundamental solutions for the free motion equation has the asymptotic form 
of the type of Hankel functions 


c£(q) = c;;r(q) ± '(Z’M ■ 


(21 


These solutions are useful for the calculation of the ^-matrix and analyzing the decaying Camov 
states. If we are interested in bound states (e = —k 2 /2,q = ik) the solution for the equation oi 
free motion with a corresponding asymptotic 

c t°r l (k) = *' \(i rta w + >('ni 9 w] ('22) 

must be used. The numerical values of solutions (16), ( 17) can lx* obtained by using the book [8], 
where the function M(a,b;z) is tabulated. Similarly to the regular and irregular solutions of the 
free motion Schroedinger equation the functions ( ' T ^ 9 and C 7 j ’ <:l are oscillating functions of n and 
the period of oscillations decreases with increasing energy e. 


3 The solution of the scattering problem in HOR 

Consider now the solution for set (12). It follows from equations (12) that the coefficient C n i for 
n > N + l obey the equation of free motion with an appropriate asymptotic, i.e. C u i = C’^, 
where C° t is the solution for the equation of free motion with asymptotic (19), (21) or (22). The 
coefficients C n i(n > N + 1 ) form the ’’external” part of the wave function in NOR. The coefficients 
C n i(n < N) belong to the ” internal” part of this function. The equation 

siintern extern 

C A II “ 1 Nt 

plays a role of 11 fitting” condition of ” internal” and ’’external” parts of the waw function. The 
r.h.s. of this equation has one of the form (19), (21) or (22). Into the left hand side of Rq. (23) 
the solution of the set (12a) must be substituted. The last one can be found in the following 
manner [2]. At first we shall diagonalize the truncated Hamiltonian matrix ||// nn '|| using the 
unitary transformation F, i.e. turn from C n i to the new coefficients 

C' xl = £; r A „C n , , A = 0, 1 , N . (24) 

n=0 

As a result of this transformation, equation (12a) takes the form 

(E\ - s)C’ xl = —VxjmTnn+iCn+u , A = 0 , \ 

! 
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i.e. 


and 


n* — 

°A / — 


£a 

Ex- e 


C 


N+ll 


N p* p 

ynr \ ' ^ An^ AN , 

Cn/ “ " ^ ~e^7' 


NN+ 1 


CV+l/ 


A=0 


where Ex - is the eigenvalues of the matrix ||//„ n '|| (n,n' < iV). 

Substituting the ” internal” solution (25b) at n = N into Eq. (23) we obtain 


N 


r* p 

n /^rO pO p \ ' An A AN 

^ NN+1 ^N+l/ “ ~ Ly Nh r nN+\ — ^ 

A=0 


E\ — £ 


Tnn 


+i 


(25a) 

(256) 


(26) 


If we deal with bound states of Gamov resonances, C® t and C% +ll are the known functions of 
energy (see (21) and (22)). In such cases the condition (26) is the transcendent equation which 
may be used to find the energies £{ of the bound or resonant states. For the scattering problem, 
we get in accordance wjth (19): 

C° Nl = C r N y + tan 6, C%r , 


C& + „ = C' N 'U, + tan 6, C™f, • 

Substituting these expressions in (26), we find in accordance with Refs. [2]: 


tan Si = — 


n Te 3 _i_ p sireg 

^Nt + rAW +! ^N+ll 

syirreg . >; s^irreg 

^ Ni • r NN+l '-'N+ll 


(27) 


(28) 


It can be seen now that the scattering phase at an arbitrary energy e can be obtained by diago- 
nalizing the Hamiltonian matrix ||// nn <|| (n,n' < N) but one time. 

In agreement with the Ritz variational principle, the negative eigenvalues Ex < 0 of the 
Hamiltonian matrix ||// nn -|| (n, n' < N) may be treated as approximate values of the energies 
of discrete levels of a particle in the studied potential. In this case the approximation accuracy 
improves with increasing the size of the' matrix N. The question arises, what is the sense of the 
matrix positive eigenvalues and of the respective wave functions? The question was answered in 
works [2, 4, 9] as follows. In the limit e — ► E\ expression (28) takes on the form 


tan St (Ex) = 


f-rreg 

i ^N+U_ 

s^irreg 

°N+1/ 


(29) 


By comparing this result with formula (19) we get the coefficient C T ^l lt (Ex) = 0 for Ex. Thus, 
by diagonalizing the Hamiltonian matrix ||// nn <|| (n,n' <• Af), we find the solutions for equations 
(12) in the region of the continuum at such discrete energies Ex > 0 which correspond to the 
vanishing of the HOR‘ wave function C n i(E\ ) at the point n — N + 1. The scattering phase 
can be calculated at such energies using simple formula (29). In the asymptotic limit of high N, 
the diagonalization of the Hamiltonian matrix on the cutoff basis n < N means the solution for 
the Schroedinger equation (9) with the additional condition Xt(b) = 0, where b = 2(N + i/) 1/2 , 
i.e. when the system is placed within a rigid box of radius 6. In this case the energy spectrum 
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for equation (9) gets discrete and the energy of any level becomes a function of the position of 
the boundary point b = 2[(N + . The condition Xi(b) = 0 is nothing other than the 

equation for the P-matrix poles in the system of radius b described by the Schroedinger equation 
(9) [10]. Thus the eigenvalues E\ of the Hamiltonian matrix ||//„ n '|| are poles of the discrete 
analogue of the P-matrix. The important point is a convergency of this approach. The practical 
calculations [3, 4] show that for rather smooth potentials V(i') it is sufficent to use a number of 
terms TV ~ 20 — 30 in the expansion (2) in order to calculate phase shifts and other scattering 
characteristics in reasonable accuracy.Therefore the using of HOR, or J-matrix approach [2], is 
a rather effective and practicable method for the study of continuum problems. Some additional 
example of application of this method will be discussed in section 6. 


4 Multichannel case 


Let us consider the case of two open (binary, spinless) channel for simplicity. The wave function 
has the form of a column 


0(r) 


M r ) 

M r ) 


(30) 


and the Hamiltonian is the operator matrix of a size 2x2: 


// 


H\\ Pl2 \ 

tin tin ) ‘ 


(31) 


Let us assume that the wave function of the entrance channel i/q(r) is characterized by the following 
asymptotic behaviour 

Mr.) ~ (e~ ikir ~ S u e lk ' T )/r (32a) 

while in the second channel only the outgoing wave presents 

Mr) ~ - ({v,lv 2 )' l2 Sne' k * r ) /r . (326) 

The transition into n-representation consists in the expansion of both channel wave functions 

V>l( r ) = EClnkroi > , 


Mr) = E Cam |m, 7-02 > 1 (33) 

m 

\ 

in terms of harmonic oscillator wave functions |n,r 0 i >, |m,r 02 > with a unique frequency hu> 
while the linear scale parameters ro; = {hj /^ca) 1 / 2 can be different for the channels 1 and 2 if the 
reduced masses fi\ and p 2 of two fragments in these channels are different. Assuming that it is 
possible to restrict ourselves to a truncated matrix of the potential energy 

Vin.in' (0 < n,Ti' < /V, ), V 2m , 2 m' (0 < m,m' < fV 2 ), 

Vln,2m , V 2 m ,1„ (0 < Tl < N x , 0 < m < N 2 ) 
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(generally speaking N X ^ N 2 ) we obtain the following set of equations for C Xn , C 2m coefficients 
instead of Eqs., (12), (23), (27) 

(H - E)C = -TC° , 

/-( f iG / 1 / '0 

< l/V, — ( IN] > ( 2/V 2 — C '2N 2 7 

Cl = Cf* - SiC+ , n > yVj , 

CL = -(qi/q^S^CL , m > N 2 . (34) 

Here C is a column of jV, + N 2 + 2 coefficients C 1 oC 1 i...C iW iC' 2 o-.C 2 n 2 , H is the matrix of 
Hamiltonian in a truncated basis |«,r 01 >, |m,r 0 2 > (n < N\, in < N 2 ). The column TC contains 
only two nonvanishing elements, namely Tn u ^i+\C^ N i+1 in the N x + 1-th row and Tn 2 n 2 +\C 2N2+1 
in the last row. The functions C% = C T ik 9 ± iC " k r * 9 are the same ones as in Eq. (21). As it 
was shown in Ref. [2, 11] the asymptotic of the function V ,± = £Lo C* |n > is of the form 
k exp(±ikr ) for — > oo. This fact and the difference of r 0t in various channels are the origin of the 
factor {q 2 /qi)^ 2 in Eq. (34) instead of the usual velocities ratio ) 1 ^ 2 in Eq. (32). Solving 

the Eq. (34) similarly to Eq. (12) we obtain the following results instead of Eq. (26) 

Cuv, = (\ni ~ Si Cj^,-, = Pu(C xNl + 1 — > c >uC XNj+1 )+ 


+Pn (— (tteA/i I'^Sif 'L 2 +i) > 

C 2N2 = -S n C+ N2 (q 2 lq x ) xli = P 2l (<SS + i - S u C+ Nl+l + 
+ P-22 (-{qi/qi ) 1/2 .s 2 ia+ V2+ ,) 


where 


P _ _ 

r *J p 


AN, 1 \Nj , 

Ex- E 


' N,N, + 1 


E\ is the eigenvalue of the truncated matrix //, (rAo-"^^---^^) * s th e corresponding eigenvector 
of this matrix. 

The relations (35) should be considered as the equations for elements of the ,S'-matrix. The 
solutions of these equations are of the form 

■Si = -£) (Ci - Ni + E nC 1Wl+1 )(C 2N2 + P 2 2C 2N2+l ) — P2\PnC 1Nl+1 C 2N2 + 1 


^ _ 1 2iP u P 2l 

‘ 21 D sfnqtf2 

D — [(CL, + P \\C xNi+x ){C 2N2 + P 2 2 C 2N2+x ) - P\ 2 P 2 \ cL, + 1 C 2N2 + 1 ] • 

Here the property of the Casorati determinant 


(36) 


Tn N+l 


c 


N 


C/v 


n- c' , + 
l N + 1 L/ n+i 


2i_ 

7 T<7 


was used. The elements S 2 2 , S X2 of the 5'- matrix can be obtained from (36) by substitution of 
indices 1 and 2. Obviously the A- matrix is symmetrical in accordance with the time reflection 
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symmetry of the Hamiltonian. As in the previous section, the eigenvalues E\ are the poles of the 
discrete analogue of the P-matrix. 

The eigenfunction^ ij)\ = £ n T** | n > are discrete analogs of "primitives” (in terms of the 
paper [10]). 

The expressions (36) allow us to find the numerical values of the A- matrix elements and then 
to calculate the cross-sections of elastic scattering and reactions 

<T' l = fi\Su - 1|'\ (37) 

= ^|S 2 ,| 2 , (38) 

the differential cross-sections, various polarization characteristics (taking into account the spin 
degree of freedom) etc. If we want to describe the reaction with three-four fragments in the final 
states it is necessary to extend the above developed formalism, which is valid only for two body 
(binary) channels, to three, four body collisions. 


5 The description of ’’true” many body scattering in a 
hyperspherical HOR 


We restrict ourselves by the case of the so called "true" many body scattering (TMBS) when the 
wave function of an A body system is in the asymptotic region of the form 


VK{P, ii) = f’K'y./voV 


-‘fcpy 


i(^) - 8' A- 'V, /Vo -vo f 

A'V 


ikpy 


r A' V 


(fi) 


OG 


(39) 


p 2 = X]/Li( r i — R) 2 is a global radius in 3 {A — 1 ) dimensional space, the angles are hyperspher- 
ical coordinates in this space. R is the center-of-mass of the system, V/,-, ( 1 1 ) is a hypersperieal 
harmonic with a global momentum A , 7 substitutes all the rest quantum numbers labelling this 
harmonic. The approximation taking into account only the contribution of TMBS is valid if then' 
is a ’’democracy” in the A body system i.e. there is no pair of particles with dominating interac- 
tion between them in comparison with the rest of the interactions. The TMBS - approximation is 
applicable to a lot of processes of three, four body decay of light nuclei and hypernuclei [17] (for 
example disintegration l2 C — » 3o etc.). 

For the description of TMBS we shall use the expansion of the A-body wave function i/>(ri ...F 4 ) 
in terms of A-l body oscillator wave function (the center of mass variable R is excluded) 

|t»A' 7 >= RnK(p)YKy(tt) (40) 

depending on hyperspherical coordinates p, f 2: 

if) = ^2 < nK^\4> > |nA'7 > . (41) 

nK 7 

Further consideration is totally parallel to sections 1-4 and we represent the result in very short 
form. Instead of Eq. (3) we have for the many body case 

Rnk(p) = p- {3A - 4)/2 <t>n(p) , ( 42 ) 
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4i{p) = (-n n < 


2 n! 


: p £+1 e-*/*Li + '' 2 (p 2 ) 


\\r(n + £ + *y 

£ = K + (3A - 6)/2, p is taken in units of r 0 . The Eq. (4) takes on the form 

y < nA' 7 |// — E\n'K'-)/' >< n K'^' \ip >= 0 . 
n'K'V 


(43) 


(44) 


The kinetic term T in the Hamiltonian H = T + V is diagonal in the quantum numbers K and 7 . 
As for the main quantum number, n, the matrix is three diagonal with respect to n and its matrix 
elements coincide with Eq. (5) except for substitution of / by £. We also truncate the matrix of 
potential energy V = Et<j V H to n < N , K < K max . Then for n > N the expansion coefficients 
< nA' 7 |V’ >= C n i obey the three-term recurrent relation similar to Eq. ( 12 b) 


sJrUn + £ + |) < n - > - (2 n + £ + -| - q 2 ^j < nK'y \xjf > + 

J{n+l){n + £+ l -)<n + 1 >- 0 , q = V2E . 

This difference equation has two fundamental solutions 


and 


CI7 = 


^ -jirreg 


2n! 


-q £+ 'e-< 2 ' 2 L £ n +1 ' 2 (q 2 ) 


l£ *C&(q) 

or the equivalent pair of solutions 

r ± _ 2( i 

U n£ — 


^r(n + £ + |) 

2 q r w)c;; £ 9 w) 

VP Jo 


q2 _ q/2 


dq' 




r 

Jo 


C^W)£2W) 

q 2 — q f 2 ± iO 


dq . 


(45) 


(46) 


(47) 


(48) 


The problem of TMBS is similar to the multichannel problem described in section 4. Thus the 
wave function ( 39 ) with ingoing wave in some channel /fo 7 o and outgoing waves in each channels 
K 1 7 ' under consideration takes on the form at n > N (in principle the truncation boundary N 
may be different in the various channels K 7 ): 


< nKn/ |V> >= , fro-ro^n /‘('Z) - S K-y,K'yC + £ ,{q) . 

A'V 


In analogy with Eq. (36) we can obtain 


S = A B 


where 


M)a'V,K7 = PK , 'y',K'yC^ +1£ (q) - 6K'y,K''y'CN £ (q) , 
= Pk'~i' ,K- y£ N + 1 /(?) — ^K-y.K'-r'C Nxitf) > 


(49) 


(50) 


(51) 


92 



^ X X\NK^_ > T 

Pl\~r,K''y' / y ¥^¥ X Nh\N+\h' 

E\ and < N K^\\ > are eigenvalues and eigenvector components of the truncated Hamoltdnian 
matrix < nK'y\H\n' I<‘ V > (n,n' < N). The poles of the .9-matrix (i.e. bound states and Gamov 
resonance states) can be found from the equation 

det .4 = 0. ( 52 ) 

Thus we have all expressions that are necessary for the construction of the wave function for few 
body states belonging to continuum or discrete spectrum in the frame of TMBS approximation. 

6 Soft dipole mode in "Li and three body continuum 

In order to demonstrate the effectiveness of the HOR in an analysis of concrete nuclear processes 
the calculations of the properties of the low energy El excitations in the 11 Li nucleus were done. 
The remarkable feature of this nucleus is a presence of a large neutron halo formed by t wo neutrons 
weakly bound with th.> 9 Li core. In 1 his connection the following model was used for the description 
of this nucleus. 


6.1 The model 

It was assumed that the "Li ground and continuum states can be interpreted in the framework 
of the three-body cluster structure 'Li+n -(- n. 

1) The cluster 9 Li is supposed to be structureless and the excitations of its internal degrees of 
freedom are not considered. 

2) We don’t account for non-central components of the interaction between two valence neu- 
trons and between valence neutron and the cluster 9 Li. Therefore, the wave function can be 
characterized by the three-body orbital angular momentum L, spin S = 3/2, total angular mo- 
mentum J and its projection M. 

3) The states with the total spin of the valence neutron pair .9 = 0 are only considered, and 
the ground state three-body orbital angular momentum is supposed to be equal to zero: L = 0. 

1) n- 9 Li interaction is described by the shallow potential of Johansen et al [13]. N N -interaction 
is described by the Gaussian potential [13]. 

5) Only democratic decay channels are allowed for. 

The wave function of the system 9 Li+n + n, y), is expanded in three-body hyperspher- 

ical functions, <S> ,X ly M (Cl) (including the internal w'ave function of J Li with a spin .9 = 3/2) 


y) = J2 'l’K ) i x iy(P)*K lvJM (P) 


where K is hypermomentum, l x and l y are the angular momenta corresponding to the Jacobi 


coordinates 




i*i + r 2 

2 


- r 3 ) , 


(54) 
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respectively, m is the neutron mass, r t are coordinates of the valence neutrons (i = 1,2) and the 
cluster 9 Li (i = 3), p— (x 2 + y 2 ) 1//2 is a three-body hyperradius. 

In the c. m. frame the Hamiltonian is of the form: 


H — T -\- \\ 2 + V13 + V23 , 


(55) 


where T is the three-body relative motion kinetic energy operator, and V tJ are the two-body 
potentials. For the radial wave functions i y (p) we have the usual set of the K -harmonic 
method coupled equations (see, e.g., [12]). The equations are solved by expanding the radial wave 
function. 


V’ft' L iy(p) = ^2 


n — 0 




(56) 


in the six-dimensional harmonic oscillator eigenfunctions. To calculate the bound state energy, 
i.e. to locate the corresponding 5-matrix pole, one should solve the nonlinear equation [1 1] 


det A = 0 , 


(57) 


where the matrix A is given by the Eq. (51). 

For the continuum spectrum states we calculate 5-matrix for any positive energy E using (50) 
The interactions of the valence neutrons with each other and with the cluster 9 Li are described 
by the potentials V u {^\ 2 ) and Fi 3 (r 13 ) = ^23(^23), respectively. We use the following parametriza- 
tion of the potentials [13]: 

VSi(r) = exp [-(r/fcjjV] + V& exp [-(r/fcg ) ) 2 ] , 
v}} ] = -31 MeV, v£ ] = 0, b\\ ] = 1.8 fm; 

= -7 MeV, U ( 3 2) = -1 MeV, = 2.4 fm, = 3.0 fm. 


In the external asymptotic region n > N we consequently allow for channels V characterized 
by K = K m i n , Kmin + 2,... ( K m i n is the minimal possible value of K for a given J) until the 
convergence for all physical properties under consideration is achieved. The convergence is found 
to be very good, and the allowance for the decay channels with K > K mxn + 2 do not yield any 
visual variation of the results. So, we consider in the external asymptotic region n > /V the 
channels with K < I\ m in + 2 only. Note, that components with all possible values of K < N are 
accounted for in the calculation of the wave function in the inner region n < N. 

The parameter hu> is set to be equal to 7.1 MeV in our calculations. This value corresponds 
approximately to the minimum of ground state energy E 0 . 


6.2 The ground state 

The results for the 11 Li ground state for different values of the truncation parameter N are pre- 
sented in the table 1. The variational ground state energies, Eq’\ obtained by the pure diagonal- 
ization of the truncated Hamiltonian matrix are listed in the second column, while the J-matrix 
results, Eq, which are the solutions of the eq. (57), are listed in the third column. It is seen, 
that by locating the ,9-matrix pole using eq. (57) that is equivalent to the allowance for the 
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Table 1: 11 Li ground state properties (see text for details). 


Truncation 
boundary N 

Ground state energy, 
MeV 

Neutron halo 
mean square radius 
< r 2 >\{ 2 , fm 

Wr~ 

E o 

< r 2 >if w 

< r 2 ■>}'’ 

12 

-0.012 

-0.150 

2.83 

3.81 

16 

-0.116 

-0.199 

2.91 

3.29 

20 

-0.171 

-0.225 

2.98 

3.31 

'24 

-0.202 

-0.240 

3.04 

3.32 

Experiment 

-0.247i0.080 

3.16T0.1 1 


long asymptotic tail of the wave function, we improve essentially the convergence for the binding 
energy. 

The results presented in the table 1 have been obtained using Lanczos smoothing of the three- 
body potential energy matrix [5, 14]. 

The "Li r.m.s. radius, < r 2 >\[\ can be calculated by the following equation: 


O CJ o '*2 

< r 2 >n= — < r 2 >9 +— <P > , 

11 11 mw 




where < r 2 ><^ 2 is the 9 Li r.m.s. radius and the mean square value of the hyperradius, < p 2 >, 
can be easily calculated using the ground state wave function. The values of < r 2 >'/ 2 ^ and 
< r 2 > 1 / 2 obtained by the pure diagonalization of the truncated Hamiltonian matrix and with the 
allowance for the asymptotic tail of the wave function, respect ively, are presented in the 1-th and 
the 5-th columns of the table 1 . 

It is seen that in calculation of the ground state, the allowance lor the wave function asymp- 
totics is very important for a weakly-bound system like "Li. 1 he terms of expansion (56) with 
the number of total oscillator quanta /V — 100 that cannot be obtained in the usual os< illatoi 
basis variational calculations, play an essential role in the description of the transverse momentum 
distribution, r.m.s. radius, etc. The convergence of < r 2 > 1/2 , transverse momentum distribution 
and other properties of the wave function in the full J-matrix calculation is rather good. Nev- 
ertheless, it is seen that the r.m.s. radius converges to a value that is somewhat larger than the 
experimental one, and the calculated transverse momentum distribution appeares narrower than 
the experimental one. These shortcomings can be overcome by the adjustment of n 9 Li potential. 
We have not aimed to fit the potential to the "Li properties, we have just take its parameters 
from ref. [13]. 


6.3 The soft dipole mode 

The dipole transition operator in our model is of the form 

M(E\p) = -^feyY l ,(y), (59) 
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where e is the proton charge, A = 11, Z = 3 and the number of valence neutrons, N v = 2. The 
operator (59) corresponds to the excitation of the three-body cluster modes only. The excitation 
energy of the first excited state of 9 Li is relatively high (~4 MeV). So, low-energy El-transitions 
correspond to the excitation of the cluster degrees of freedom only and should be described by the 
operator (59). 



E, MeV 


Figure 1: Comparison of our results for B(E\\g.s. — > continuum) in 11 Li with results of other 
authors. 1 this work (J - matrix method), 2 — ref. [15]., 3 ref. [16], 4 -- experimental data 
paramet ri/al ion of ref. [17]. 

The cluster reduced probability of the El -transition, B(E 1; Ej — E 0 ), associated with the 
operator (59), is displayed on the figure 1. 

This figure shows the comparison of the results of our calculations of cluster B(E 1; Ej — Eo) 
with the parametrization of experimental data of ref. [17]. The agreement is reasonable. The 
form of the B(E 1; Ej — E 0 ) peak is well reproduced, the discrepancy in the position of the 
B(E 1; Ej — E 0 ) maximum is supposed to be eliminated by the adjustment of the potentials. The 
results of the B(E 1; Ej — Eq) calculations of refs. [15, 16] are also depicted. All these colculations 
give a low energy peak which can be associated with the soft dipole mode. 

The soft dipole mode exhausts about 90% of the cluster sum rule (EWSR) associated with the 
operator (59). The contribution from the soft dipole mode to the total EWSR is relatively small. 
In the vicinity of the sharp B(E 1; Ej — E 0 ) maximum at the excitation energy E ~l-2 MeV only 
~8% fraction of the total EWSR is exhausted. Nevertheless, the account for the soft dipole mode 
results in an essential increase of the electromagnetic dissociation cross section of 0.8 GeV/nucleon 
"Li beams on Pb and Cu targets. Using the sums of the "Li and target nucleus charge radii as 
impact parameter we obtain for the electromagnetic dissociation cross sections the values of 0.966 
barn for the Pb target and 0.132 barn for the Cu target; the corresponding experimental values 
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are 0.890 ± 0.110 batn and 0.21T0.04 barn, respectively [18]. E 0- and E2-transitions give only 
1.2% contribution in the cross sections. 

Thuse, it is shown, that cluster model 9 Li+n+n yields a good description of the ground state 
properties and El-transitions in the 11 Li nucleus. The HOR may be used successfully in the 
studies of weakly-bound systems with long-tailed wave functions, e.g., in the study of neutron 
halo properties. For both bound and continuum states the correct account of the wave function 
asymptotics in the framework of the oscillator representation of scattering theory is very important 
in such studies. Low-energy El-transitions in "Li are of the cluster nature. The widths and the 
position of resonant states calculated in the democratic decay approximation are in a reasonable 
agreement with experiment. 

Appendix. Isolated States. The scattering problem with nonlocal separable potential V A 
can display some peculiarities which we explain here using a simple example when t he Hamiltonian 
H is approximated by the matrix of a size 2x2 (i.e. N = 1 ). In specific situation when / ,>i = — Voi , 
i.e. the nondiagonal matrix elements of the kinetic and potential energies cancellate each other 
Ho, = 0, we obtain that the Harmonic oscillator wave function /?oo( r ) an eigeniumt ion of this 
Hamiltonian corresponding to the eigenvalue Eq = T 0 o + loo- If E 0 > 0 we find an example of the 
bound state embedded in continuum [19]. It is clear that the eigenfunction R m (r) is not connected 
with the rest basis states R n o(r). Thus it is isolated from continuum states and can be calk'd an 
isolated states. The phase shift 6 0 (k) displays a narrow resonance near energy E\, at. small value 
of ho, ■ It transforms into the resonance of zero width when Ho, — > 0. 
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